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Abstract 

We in this paper study theoretically the spin-polarized quantum transport through a T-shape 
quantum dot-array by means of transfer-matrix method along with the Green's function technique. 
Multi-magnetic fields are used to produce the spin-polarized transmission probabilities and there- 
fore the spin currents, which are shown to be tunable in a wide range by adjusting the energy, 
and the direction-angle of magnetic fields as well. Particularly the opposite- spin- polarization 
currents separately flowing out to two electrodes can be generated and thus the system acts as a 
spin splitter. 
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INTRODUCTION 



Recent progress in tlie nano-fabrication of quantum devices enables us to study electron 
transport through quantum dots(QDs) in a controllable way since a wide range 

of energy can be achieved here by continuously changing the applied external potential in 
contrast with real atoms. A QD- array regarded as an artificial crystal^, ^,0] also becomes 
reliable due to the recent development of nanotechnology. When the size of structures is 
small with respect to the coherence length, the transmission probability plays an essential 
role in the quantum transport. For example, we can obtain the conductance of the devices 
from the transmission probability with the help of the Laudauer-Biiticker formula. Various 
Deen developed for numerical calculations of the transmission probability 



methods have 



system [23, 
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3, 3: 



17 



18|. The quantum transport through a multi-terminal 



25| can be analyzed in terms of all transmission probabilities between any pair 



of terminals in principle. Biiticker et al.[19| derived a scattering matrix for a three-terminal 
junction consisting of homogeneous wires based on the unitary condition and the assumption 
that the scattering matrix is real. It was shown that this matrix is energy independent and 
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the coupling parameter has to be determined from phenomenological arguments 20|, |2l|. An 
energy-dependent scattering matrix for three-terminal junction was found later [22 1. 

Spintronics has become a new branch of condensed matter physics and material science 
where not only the charge but also the spin degree of freedom of electrons play an essential 
role26j. It is of fundamental importance to study the mechanism to generate and control 



the spin-polarized currents in spintronic devices which may promise practical applications 



in quantum computing and information 



28( 1 ■ In recent yeas the spin-polarizec 



in quantum dots have attracted considerable attentions both theoretically 
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32,133, 



and experimentally|4ll. l42l] . where the spin currents are generated 
either by magnetic material electrodes or by the rotating-magnetic field which results in 
the spin-flip. Most of the studies are concentrated on the one or two-QD devices with only 
one output electrode and the spin-polarized transport in QD-array consisting of arbitrary 
number of QDs with multi-terminal has not yet be investigated. In this paper we study 
the quantum transport through a T-shape QD-array with three electrodes in which the 
spin-polarized current is induced by static magnetic fields applied on the three arms of the 
T-shape QD-array. We calculate the transmission probability in terms of the scattering 
matrix which is spin-polarization dependent in the presence of magnetic fields and leads to 
the spin-polarized currents. It is the main goal of the present study to show how the spin 
currents in output terminals can be generated and controlled. Particularly two spin currents 
of opposite-spin-polarization flowing simultaneously out to electrodes are obtained and seem 
extremely useful in the design of quantum- logic-gate devices. 



II. MODEL AND FORMULATION 



The system consisting of three perfect wires (electrodes) attached to a T-shape QD-array 
of arbitrary number of QDs is explained schematically in Fig.l. The horizontal QD-array is 
denoted by the lattice sites (-A^l, 0), (-A^^ + 1, 0), (0, 0), {Nr -1,0), {Nr, 0) and the 
vertical QD-array by the lattice sites (0, 1), (0, 2), (0, Njj) where A^^, Nr, and Nu are the 
numbers of QDs on the left, right and vertical arms respectively. Three perfect wires are 
also described in terms of discrete lattice sites. The first magnetic field in the z-direction is 
assumed to be applied ,for example, on site (— A/'l,0), the second magnetic field on the site 
{Nr^ 0) is along a direction of angle 6i with the z-axis, while the third one on the site (0, Nu) 
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has a direction angle 62. The latter two magnetic fields not in the same direction as the first 
one result in the spin-fiip of electrons, which is the mechanism generating the spin-polarized 
transport in the QD-array (see the Hamiltonian below). We remark that it is not necessary 
to have the three fields to apply on the certain single QD and what we choose the three 
QDs is only for the convenience of analysis. It may be more easy in practical experiments to 
superpose the external fields on the three arms of the T-shape QD-array. The Hamiltonian 
of the system includes five parts 

where 

Nr-1 
n=-NL+l,a 

+ (^iVn + <79/^H -Si cos6'i)(i+ ^ d,„ + qa„ BisinOiid'^ ^ + h.c) 

Nr-1 

— > (tdt , d, + h.c) — (t d'j' , d,., „, -|- h.c) — it d'f , d, -|- h.c) 

/ J V (n+l,0),o- (n,0),cr / V (jVjj-l,o),CT (JV^,0),CT / V (_Ar^_|.i^o),CT (-Ni^,0),a > 

n=-NL+l,a 



and 



Nu-l 

Hd,v = ^ ^0,„^^,„),^<^(0,n),a 

n=0,a 

+ (£o,Nu + (7gi^B,B2 COS ^2)c^J_^^,_/(o.^„).. + g^^B^B2 sin ^2(c?+ ^^)_/(o,^„),_. + h.c) 



n=0,(T 



itd'^ d,, , + h.c) - (td+ d,, „ 1^ + h.c) 



denote the Hamiltonians of the horizontal and vertical QD-arrays respectively with e„ and 
£o,n being the energy eigenvalues of corresponding QDs. c/^ ^^^^ '^^ „) ^ (with a =|, J,) are 
the creation operators of the electrons with spin index a on the QDs, where —a denotes the 
opposite spin polarization with respect to a. The Zeeman terms induced by the external 
fields on the sites {N^-,^), and {Q,Nu) result in the spin-flip. It should be emphasized 
that the 9i and 62 dependent terms describing the spin-flip in the Hamiltonian are the key 
mechanism of the spin-polarized transport in our system. The matrix elements defined by 
(n -I- 1, 0| |n, 0) = (0, n -I- 1| iy |0, n) — t denote the hopping integrals between the nearest 
neighbors of QDs and are independent of spin polarization , except for n — Nr , n — —N^ — l, 
and n — Nu, where Zeeman energy induced by the external magnetic fields leads to the 



energy-level splitting and therefore the imbalance of tunnel couplings between spin-up and 
spin-down electrons. The hopping matrix elements in connection with the QDs {Nji,0), 
{-Nl, 0) and (0, Nu) are assumed to be {Nr, 0| H {Nr -1,0) = (-A^^, 0| H \~Nl + 1,0) = 
(0, Nul H \0, Nu — 1) = t' . The Hamiltonians of horizontal and vertical electrodes are 
written as 

Heih— > SfO,^ , a, — > (ta^ a, „, + h.c) 

ei,n / J L (n,0),<T (n,0),o- / J V (n+l,0),o- (»i,0),<T / 

n<-(iVi+l),(T n<-(Ari+2),(T 

+ > £„a^ a, n^ — > {t<if , a, + h.c), 

/ J R (n,0),a (1,0), <T / > V (n+l,0),<T ('»,0),<t 

n>JVH+l,o- n>N 11+2,(7 

Heiv= > £rrO^ o,„ , — > (ta^ , a,„ , +/i.c), 

e',^' / ^ t/ (0,n),a (0,n),CT / > ^ (0,n+l),a (0,n),<T /' 

n>JV(/+l,CT n>Nu+2,a 

where £^ , £^ and are the on-site energies in the right-hand-side, the left-hand-side and 
the vertical electrodes respectively, ^ and ^ are the creation operators of electrons 
in the electrodes. Finally the hopping of electrons between the QD and electrode is described 

by 

where Vr^^, Vl^o-, Vu,a respectively denote the tunnel couplings between the three special 
QDs, where the magnetic fields are applied, and corresponding electrodes. 

The main quantities which we have to calculate in our formulation are the transmission 
probabilities Tij of electrons from electrode-j to electrode-z, and the reflection probabilities 
Rjj in the elcctrodc-j which is considered as input terminal. For our three-terminal system 
we have j — 1 (the left electrode) and i — 2,3 (the right and vertical electrodes) with T21 
denoting the transmission probabihty from the left to right electrodes and T31 from the left 
to vertical electrodes. However the situation is not simple for the evaluation of transmission 
and reflection probabilities in the multi-terminal system since the electron transport between 
two electrodes is by no means an isolated process but affected by other electrodes. To this 
end we begin with the stationary Shrodinger equation 

if !*) = £; I*), (2) 

which possesses a general solution of the form 

00 00 
I*) = E ^(n.o),.<,o),. |0) + E C(o,n),.a+ |0) . (3) 



Inserting Eq.(3) into the Shrodinger equation (2) yields 

(E — en,o)C(n,0),a + tC(^n-l,0),a + tC(^n+l,0),a = 0, n 7^ —Nl, 0, Nr, 
{E - £o,n)C'(0,n),cr + iC'(0,n-l),(T + ^C'(0,n+1),(T = 0, 77. 7^ Nu , {n > 1) 

{E - £o,o)C'(o,o),<T + ^C'(o,i),<T + ^C'(-i,o),<T + ^C'(i,o),o- = 0, n = (4) 

-Nl-1,o- + t'C-NL+l,<7 — 0, (5) 
{E - (ffATfl.O + CrgfJ'BiBi cos Oi))C(^NR,0),a + Bi sin OiC(^NR,0),-a 



(6) 



(7) 



{E - {eo,Nu + (TgfMB2B2 COS 92j)C(o^Nu),a + ^2 sin e2C(o,Nu),-a 

Notice that in the input electrode- 1 the wave function is a superposition of incoming plane 
wave of unit amplitude and a reflection wave with amplitude rii((j), while in the output 
electrodes-2, 3 the outgoing plane waves possess transmission amplitudes t2i(c) and t3i(cr) 
respectively, we have 

C(n,o),. = e^^""« + rn((T)e-^^>", n < -{N^ + 1), (10) 

C(„,o),. = t2i(f7)e^^>", n>NR + l, (11) 

Cio,n),a = hi{a)e'''^^\ n>Nu + l, (12) 

where, K^, K^, and are the wave vectors in the left-, right-, and vertical-electrode 
respectively. Using these wave functions, we can eliminate all the coefficients, {C{Q^n),a}-i in 
Eg. (4) and Eg. (7) and obtain 

{E - enfi)C{nfi),a + tC{n-l,Q),(T + *C'(n+l,0),CT = 0, n 7^ --/Vl, 0, Nr, 

{E — (£0,0 + So!o)^(o,o),(T + ^C'(o,i),(T + ^(^(-i.o),^ + ^C'(i,o),cT = 0, n = 

{E - (£-JVi,0 + Crfl'/^S-S))C(_Ar^,o),<T + Vi„aC-JVi,-l,CT + t'C_NL+l,a = 0, 

(^ - (sat^.o + crg/iBiBi cos 6'i))C(iv«,o),o- + -Bi sin 6'iC(iv«,o) 



-Ft'C(jVH-l,0),(T + VR^aC{NR+l,0),a — 0, 



(13) 
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where £0,0 + £0,0 effective energy on site (0, 0) witli 



<o = 7^- ^^^^ (14) 



^—£0,2- 



it'? 

E-eo,Nu-T.'' {E)- 

Here ^ {E) = y^^G'^(£;) denotes tfie self-energy due to the couphng with the horizontal 
QD-array, and G^{E) is the Green function of the lattice-site (0, Nu) satisfying the recursive 
equation 

G''{E) = [E-eu-t'G''{E)]-\ (15) 
the solution of which is seen to be 



G'^iE) = —{{E - eu) - - {E - sufY^}. 



(16) 



We can see that, the Eq. (13) does not include any of the expansion coefficients, {C(o,n),(T}, 
associated with the vertical lattice-site. Thus it imphes that the three-terminal device 
reduces to an effective two-terminal system and the effect of the vertical QD-array appears 

Eu 
{E) which has now been included in the effective on-site energy. 

Now we rewrite Eq. (13) as a matrix equation 



C(n+l,0),cr 

C(n,0),(T 
C(n+l,0),-o- 
C(n,0),-a J 



M[(n,0),E] 



C(n,0),o- 

C(n-l,0),cr 
C(n,0),-(T 
C(n-1,0) -a J 



where M[{n,Qi),E\ is called the transfer matrix which hnks the expan- 
sion coefficient vector (C(„+i,o),<7, C'(n,o),a, (^(n+i.o) -a, (^(n.o) -a)'^ to the vector 
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(C(n,o),(T, C^ri-i,Q),a, C(n,o) -a, C(„_i,o) -a)^ and thus is defined by 



M[{n,0),E] 



M[{Nr,0),E] 



t 


-1 












1 





















E-e 


n,0 2 







t 












1 


_ 








,0 


t' 


g/LiSj Bi sin 0\ 









Vr,„ 


VR,a 




1 















gUBiBi sin 6^1 







t' 










VR,a 










1 








where e'^^ is the effective on-site energy given by 



(17) 



4,0 = ^n,o, n ^ -Nl, 0, A^R 

4,0 = ^0,0 + ^ = 0) 
^'iVfl,o = ^iVfl,o + crg/J'S^Bi cos 9i, n = -Nr. 



(18) 



From Eqs. (10) and (11), it can be shown that the transmission amplitude t2i{cr) is related 
to the incident coefficients via the equation 



" t2l{<j) ' 




1 





= T{E) 






t2l{-(T) 




1 







.rii(-(7) _ 



(19) 
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where the transfer matrix T(E) for the whole system is given by 

' ^-iKR{Na+l)a 



T{E) 



e^^-^'^^^+i)" 













-iK^^{Nii+l)a 





JKR^{Nii+l)a 



X 













1 


1 
























1 


1 



-{Nl+1) 

n M[{n,0),E] 

n=Nii+l 



X 



1 





















-iJs:^(Ari+i)a 











,ii^^(Afz,+l)a 









^-iKt^{NL+l)a 





Thus from Eq.(19) we can obtain the reflection amplitudes and rii(J,) as 

rii(T) = 



(Til _l Zki^ -l ( t^-l , T43 'i 
T22 ^ T22' ^ T22 ^ T44 T44 / 



rii(i) 



-( 



T22 T44 



T4I I T43 '\ _i_ T42 ^ T21 I T23 



_|_ £43 'J _j_ -'42 t 
T44 T44 ' T44 V T22 



1 



^24 ^42 
222 T44 



and the transmission amplitudes t2i(T) ^iid t2i(i) 

t2i(T) = Tn + Tisrnd) + Tis + 7\4rii(i), 

t2i(i) = Tai + r32rii(t) + Taa + r34rii(i). 
The spin-polarized reflection and transmission probabilities are seen to be 

and the total transmission and reflection probabilities are defined by 







(20) 



(21) 
(22) 

(23) 
(24) 

(25) 
(26) 



T21,tot=T2i(T)+T2i(i), 



(27) 
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Rll,tot^R2l{'i)+R2l{i). (28) 

The spin current can be derived from tiie difference between tiie two spin-polarized trans- 
mission probabilities 

AT2i = T2i(T)-r2i(i). (29) 

The transmission probabihty, T31 (cr), from electrode-1 to electrode-3 can be simply obtained 
by replacing the on-site energy s'q^q = £o,o+£o!o matrix M[(0, 0), E] with s'q q = £0,0+^0^0 



where 

- 

Cfi n — 

E — £1,0 



4o - (30) 

E-£2,0 



it'? 

and (E) = Vl^G^{E) is the self-energy due to the presence of the right electrode, and 

G'^iE) ^[E-eu- t^G'^{E)\-\ (31) 

is the Green's function of the lattice-site {Nji -|- 1,0). Moreover the effective on-site energy 
£^ Q of Eq. (18) should be replaced by 

4,n = ^o,n, -Nl,0,Nu 
^'-Nl,o = ^-Nl,o + ctqubB, n = -Nl 

^0,0 = ^0,0 + ^0^0) = 0, 
^0,Nu = ^0,Nu + (^9fJ'B2B2 COS 02, U ^ Nu ■ 

Then the transmission probability r3i(cr) can be evaluated from 

(32) 

and finally we obtain the total and spin transmission probabilities T^i^toti ^T^i in the same 
way as for the calculation of T2i,tot and AT21. 

The conductance of the device can be obtained from the transmission probability with 
the help of the Laudauer-Biiticker formula 

G = (^)T, (33) 
10 



where e is the electron charge, h the Planck's constant, and T the transmission probability 
of the device. 

III. NUMERICAL RESULT AND DISCUSSION 
A. Homogeneous dot-array 

In the numerical evaluation the lattice constant a is chosen as unit of length and hopping 
constant t as the unit of energy. The on-site energies of all QDs, as well as Si, Sr , Su are set 
to the same value of 2t. The coupling parameters between the QDs and the three electrodes 
are assumed to be Vl,(j — Vr^^ — Vu,a — 0.5t and q^jle^Bi — g^B2B2 — li- 

First of all we consider the simplest case that there is only one-QD on each arm of the T- 
shape array besides the central dot, i.e. = Nr = Njj = 1. Figs. 2(a), 3(a) and 2(b), 3(6) 
display the energy-dependence (spectrum) of the spin-up and the spin-down transmission 
probabilities to the electrode-2 and the electrode -3 respectively (for t = 1 , t' = 0.8, 6i = jTT, 
02 — jn). The total transmissions T2i^tot , 2^3i,tot are shown in Figs. 4(a) and 4(6). It can be 
seen that in the presence of the external magnetic fields, the transmitted electrons split into 
spin-up and spin-down components with different spectra. Fig. 5 is the plot of the relative 
transmission probabilities between two spin components AT21 ( solid line) and AT31 (dotted 
line) versus the energy respectively. We can see that in some electron-energy ranges, for 
example from 1.8t to 2At or from 3.4t to 3.8t, the spin transmission probability to the 
electrode-3 ,Ar3i, is positive (i.e. the spin current with spin-up polarization) while the spin 
transmission probability to electrode-2, AT21, is negative (namely the spin-down current). 
However the situation would be opposite in other energy region , for example from Ot to 1.2t 
, where AT'21 is positive and AT31 is negative. 

Figs. 6 show the spin transition probabilities AT21 (solid line) and AT31 (dotted line) as a 
function of energy for the multi-QD case that A^^ = Nr = 4 and Njj = 3 with the magnetic 
field direction angles 61 — 62 — In the electron-energy range from 0.6t to 0.75t or from 
l.T&t to 1.90t the spin transmission probability ATsi is positive (i.e. the spin-up current 
flowing to electrode-3) while AT21 is negative ( spin-down current flowing to electrode-2). 
Opposite spin currents are observed in the energy region from 0.38/; to 0.47/; , where ATsi is 
positive and AT21 is negative. Fig. 7 are the results of AT21 (solid line) and AT31 (dotted 
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line) for the direction angles 9i = jTt while 62 = |vr. We see that the spin currents are 
sensitive to the fields. It may also be worthy to point out that in the energy region between 
2.9t and 3. It, AT21 is negligibly small while AT31 is positive. In other words the spin current 
can exist only in one output electrode. Comparing with previous case (one-QD on each arm) 
more resonant peaks of the spectra appear with the increasing number of QDs in the same 
way as the usual quantum transport in the QD-array without spin-polarization. 

The generation of spin currents in the two output electrodes in our model can be under- 
stood that a force applied on the spin-S which is induced by the spatially inhomogeneous 
magnetic fields, namely F = g'/ibS ■ yB ( where fib is the Bohr magneton, g is the spin 
g- factor), removes the degeneracy of spin-polarization in the transport and thus leads to the 
relative shifts of resonant peaks of the transmission spectra between the spin-up and spin- 
down components of electron. The spin current appears when the transmission spectrum of 
one spin component reaches a maximum while the spectrum of the other spin component is 
a minimum. 



B. Effect of disorder 

The homogeneous dot-arry is just a theoretical idealization and disorder is necessary to 
be taken into account in practical systems. It is interesting to see the disorder effect of the 



QD-array on t 
following Refs. 



le spin-po 



43 



44 



45 



arized 
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quantum transport. To this end we consider the disorder 
47( 1 that the on-site energies of QDs are alternated with 4(1 — x) 
and 4 ■ x, where < x < 1 is the probability distribution parameter. When x = 0.5, the 
system reduces to the non-disordering one. In the T-shape QD-array on-site energies ^ 
{—Nl < n < Nji , on the horizontal arm) and Eg^ {I < n < Nu , on vertical arm) are set 
to 4(1 — x) for odd-Tj, and 4x for even-Ti. Fig. (8) is a comparison of plots of AT31 for the 
probability distribution parameter x = 0.47 (dotted line), 0.5 (solid line) respectively with 
A'^^: = = 4 and A*"^ = 3. We can see that the disorder can shift the resonant peaks of the 
transmission probability AT31 slightly, and suppress the height of some peaks but no more 
then 20% while increase the height of other peaks no more then 18%. In other words the 
spin splitting is robust against the disorder. ^ 

The latest advances in nanotechnology make it possible to fabricate QD-array[3] and 
the model of spin-splitter proposed in this paper may be realizable experimentally. In 
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practical experiment it is more easy to apply the external fields on three arms of the T-shap 
QD-array instead of three QDs. The relative shifts of resonant peaks of the transmission 
spectra between the spin-up and spin-down components of electron can be modulated by 
the magnetic fields and on-site energies. 

IV. SUMMARY AND CONCLUSION 

In this paper we demonstrate theoretically a method to generate and control spin cur- 
rents in a three-terminal system. The output spin currents are tunable in a wide range of 
magnitudes and various output configurations by adjusting the energy and the direction an- 
gles of the magnetic fields as well. Moreover it is demonstrated by the numerical evaluation 
that the spin currents remain in the presence of the disorder of QD-array. Particularly the 
spin currents with opposite spin polarizations in the two output electrodes can be gener- 
ated. Thus this device is, as a matter of fact, a spin splitter similar with the light beam 
polarimeter in optics. This observation may have practical application in the fabrication of 
spintronic devices for logic gates with spin-up and spin-down outputs regarding as qubits. 
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FIGURE CAPTIONS 

Fig.(l) T-shape QD-array of three-electrode with three magnetic fields applied on the 
dots {-NL..O),iNR,0) and (0, A^c/) respectively. 

Fig. (2) Spin-up (a) and spin-down transmission probabilities to the electrode-2 as a func- 
tion of energy. 

Fig. (3) Spin-up (a) and spin-down ransmission probabilities to the electrode-3 as a func- 
tion of energy. 

Fig. (4) Total transmission probabilities to electrode -2 and electrode- 3. 

Fig. (5) Spin-polarized transmission probabilities to electrode- 2 and electrode-3. 

Fig. (6) Spin-polarized transmission probabilities to electrode-2 (solid line) and electrode-3 
(dotted line) versus the energy for angles 9i = ^tt and 62 = |vr. 

Fig. (7) Spin-polarized transmission probabilities to electrode-2 (solid line) and electrode-3 
(dotted line) versus the energy for angles 9i = |7r and 62 = |vr. 
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Fig. (8) The comparision of the plots AT31 with (dotted hne) and without (sohd hne) the 
disoder for angles 9i — ^tt and 62 — ^tt. 
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